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Abstract: The method of analytic continuation from imaginary to real chemical potential 
is one of the most powerful tools to circumvent the sign problem in lattice QCD. Here we 
test this method in a theory, two-color QCD, which is free from the sign problem. We find 
that the method gives reliable results, within appropriate ranges of the chemical potential, 
and that a considerable improvement can be achieved if suitable functions are used to 
interpolate data with imaginary chemical potential. 
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1. Introduction 

The phase diagram of QCD in the temperature - chemical potential plane is the subject of 
many present investigations. Understanding the different phases and the transitions among 
them has strong implications in cosmology, in astrophysics and in the phenomenology of 
heavy ion collisions. Unfortunately, perturbation theory and approaches based on effective 
models can handle a limited number of issues of the QCD phase diagram and, in fact, the 
lattice formulation is the only tool for a quantitative approach to the problem based on 
first principles. For non-zero chemical potential, however, the QCD fermion determinant 
becomes complex and the probability interpretation of the QCD Euclidean action, neces- 
sary for the standard Monte Carlo importance sampling, is lost, this being the well-known 
"sign problem". 

Several methods have been invented to circumvent this problem (for a review, see [1] 
and [2] ) : the reweighting from the ensemble at /i = [3-6] , the Taylor expansion method [7- 
15], the canonical approach [16-19], the density of states method [20-26] and the method 
of analytic continuation from an imaginary chemical potential [27-48]. Their application 
has allowed to get relevant information on the critical line separating the hadronic phase 
from the quark-gluon plasma phase in the region fj,/T < 1. 

In this paper we focus our attention on the method of analytic continuation. The 
idea behind this method is very simple: numerical simulations are performed at imaginary 
chemical potential, [i = ifij, for which the fermion determinant is real, then Monte Carlo 
determinations are interpolated by a suitable function and finally this function is analyt- 
ically continued to real values of [i. This method is rather powerful since the coupling (3 
and the chemical potential fi can be varied independently and there is no limitation from 
increasing lattice size, as happens with other methods, like those based on reweighting. 
There is, however, an important drawback: the periodicity of the QCD partition function 
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Figure 1: (Left) Phase diagram in the (T, 8) plane according to Ref. [49]. (Right) Tentative phase 
diagram in the (T, 0) plane after the inclusion of the chiral critical lines. 

and the presence of non-analyticities arising for imaginary values of the chemical poten- 
tial [49] restrict the region useful for numerical determinations to the strip < [ii/T < tt/3, 
or even less in presence of "physical" phase transitions. This implies that the accuracy in 
the interpolation of the results at imaginary chemical potential has a strong impact on 
the extension of the domain of real fi values reachable after analytic continuation. In that 
sense it is very important to answer the question about which is the optimal way to extract 
information from data at imaginary chemical potential, i.e. which is the best choice for 
the interpolating function, which only in a some cases can be guided by physical intuition, 
leading to some particular prediction for the behaviour at real [i. Moreover one should 
always be careful about the actual ranges of applicability of the method, which can be 
influenced by the various physical and unphysical transitions present in the QCD phase 
diagram, leading to possible non-analyticities. 

So far, the method of analytic continuation has been applied in SU(3) with nj = 
2 [29,30], nj = 3 [31,32] and rif = 4 [33-35,43]. Moreover, it has been tested in several 
theories which do not suffer from the sign problem, by direct comparison of the analytic 
continuation with Monte Carlo results obtained at real fi [28,36,37,44]. In most of these 
applications, a truncated Taylor series (or, more simply, a polynomial) has been used as 
interpolating function, sometimes a Fourier sum for the low temperature region [33,34]. 

The aim of this paper is to study limitations and possible improvements of the method 
of analytic continuation, by considering its application to SU(2) or two-color QCD. This 
theory is free from the sign problem and Monte Carlo numerical simulations at real values 
of the chemical potential are feasible. This allows to compare the extrapolations from 
imaginary to real chemical potential with direct determinations allowing at the same time 
both to discriminate among different Ansatze for the interpolating functions and to directly 
test the range of reliability of the method itself. The experience gained in this way can 
then be hopefully used as a guide in applications to the real theory. 

The paper is organized as follows: in Section ||, we briefly recall some general properties 
of the phase diagram of SU(N) gauge theories in the temperature - imaginary chemical 
potential plane and discuss their implications on the method of analytic continuation; 
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Figure 2: Phase diagram in the (/?, /2/)-plane; AT is the number of colors, N T the extension of the 
lattice in the temporal direction. The numerical values for (3e and (3 C are valid for SU(2) in presence 
of n/ = 8 degenerate staggered fermions with mass am — 0.07. 

in Section ^, we present our numerical results and discuss both the choice of the best 
interpolating function, showing that functions different from polynomials can considerably 
improve the method, and the ranges where analytic continuation is reliable; finally, in 
Section f|, we draw our conclusions. 

2. Theoretical background 

Long ago Roberge and Weiss (RW) have shown [49] that the partition function of any 
SU(N) gauge theory with non-zero temperature and imaginary chemical potential, fi = ifij, 
is periodic in = m/T with period 2n/N and that the free energy F is a regular function 
of 9 for T < T E , while it is discontinuous at = 2n(k + l/2)/N, k = 0, 1, 2, . . ., for T > T E , 
where Te is a characteristic temperature, depending on the theory. The resulting phase 
diagram in the (T, #)-plane is given in Fig. |l| (left), where the vertical lines represent first 
order transition lines. This structure is compatible with the /i — > —fi symmetry, related 
with CP invariance, and with the Roberge- Weiss periodicity. The ///-dependence of any 
observable is completely determined if this observable is known in the strip < < tt/N. 
It may be useful to recall the two steps in the proof of periodicity in SU(N): first, the phase 
transformation 

ip(x,r) — > exp(irfj,i)ip{x,r) , (2.1) 

then a gauge transformation with periodicity up to an element of the center group Z(N), 
i.e. a transformation with gauge group elements U(x, r) satisfying the boundary condition 

U(x, aN T ) = exp(2mk/N)U(x, 0) , k integer , (2.2) 

where N T is the lattice size in the temporal direction and a is the lattice spacing. The RW 
periodicity of the partition function extends to the observables which are left unchanged 
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by these two transformations. This is certainly the case of the chiral condensate {ipip)- 
The Polyakov loop L = Yl T= i n t U (x , ar) , instead, takes the factor exp(2irik/N) under 
the transformation fl2.2D, which implies that (L) moves, continuously or discontinuously 
according to the temperature, from one Z(N) sector to the other when (ij passes from one 
RW sector to the next. As a consequence, the chiral condensate has the same periodicity 
in ni/T of the partition function, 2ir/N, while (L) has periodicity in ni/T equal to 2tt. 
These predictions have been confirmed numerically in several cases [29,30,33-37]. 

A phase diagram like that in Fig. |l] (left) would imply the absence of any transition 
along the T axis in the physical regime of zero chemical potential for any value of N, 
nf and the quark masses, which cannot be true. Therefore, it is necessary to admit that 
the phase diagram in the (T, 0)-plane is more complicated than in Fig. |l| (left). The 
simplest possibility is given in Fig. [I] (right), where the added lines generally represent 
transitions which can be first order, second order or crossover, and can be considered as 
the continuation of the physical critical line taking place for real chemical potentials. The 
temperature T c is the critical or pseudo-critical one for the transition at zero chemical 
potential. The temperature represents the endpoint of the RW transition lines: the 
fact that the continuation of the physical critical line ends right on T E is not expected a 
priori, but is the result of numerical investigations [29,30,33-35]. 

It is convenient to redraw the phase diagram of Fig. [l] (right) in the (/3, /t/)-plane 
(Fig. H|), where (3 = 2N/g 2 , fij = afj,j is the imaginary chemical potential in lattice units 
and it has been used the fact that T = l/(aN T ). 

Given this phase diagram, it is possible to distinguish three different regimes, corre- 
sponding to different ranges of temperature (i.e. of /?), where analytic continuation can 
apply differently. 

Regime a: T > T E (or (3 > E ). 

This regime corresponds to temperatures for which the only expected non-analyticity 
at imaginary chemical potential is represented by the RW transition line. In this case the 
useful interval in fij for numerical simulations is [0, vr/8]. On the side of the real chemical 
potential, no transition line is expected. This situation is, in some sense, the best possible 
for the application of the method of analytic continuation. Simulations at imaginary \x can 
be done on a relatively large interval and, if the optimal interpolating function is found, 
its continuation should reproduce data for any real value of [i. The last expectation could 
actually be wrong if the critical behaviour induced by the RW line had some influence 
on the ranges of analyticity for the partition function also for real values of the chemical 
potential: this is an important point that can be directly checked in two-color QCD. 

Regime b: T c < T < T E (or C < p < (3 E ). 

This regime corresponds to temperatures for which a non-analyticity is expected at a 
\±i value smaller than ir/8. On the side of the real chemical potential, no transition line 
is expected. This situation is similar to the previous, with the important difference that 
the useful interval in fij for numerical simulations is restricted and the critical behaviour 
induced by the transition line may be different, thus making in practice more difficult to 
find the optimal interpolation. 

Regime c: T < T c (or p < (3 C ). 
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Figure 3: Negative side of the horizontal axis: imaginary part of the fermionic number density vs. 
the imaginary chemical potential at (3 = 1.90. Positive side of the horizontal axis: real part of the 
fermionic number density vs. the real chemical potential at (3 = 1.90. The green (blue) solid line 
represents the polynomial (ratio of polynomials) interpolating function; the dashed lines give the 
corresponding uncertainty, coming from the errors in the parameters of the fit. 

This regime corresponds to temperatures for which no non-analyticities should be met 
in fij. This implies that /}/ can be varied at will, although no additional information for the 
observables of interest here can be gotten by going farther than /}/ = 7r/4, owing to the RW 
periodicity. This regime of temperatures is probably the most interesting for physics, since 
a transition is expected here for a certain real value of the chemical potential. This implies 
that, no matter how good is the interpolation of data at imaginary chemical potential, its 
continuation to real /i should fail to reproduce data above a certain value. 

3. Numerical results 

We have performed numerical simulations on a 16 3 x 4 lattice of the SU(2) gauge theory with 
nf = 8 degenerate staggered fermions having mass am = 0.07. For this theory the tentative 
phase diagram looks like in Fig. [2], with (3e — 1-55 [36,37] and (3 C ~ 1.41 [50]. The algorithm 
adopted has been the usual exact 4> algorithm described in Ref. [51], properly modified for 
the inclusion of a finite chemical potential by multiplying the forward (backward) temporal 
part of the Dirac matrix by (e -/i ), for the case of a real chemical potential, and by e 1 ^ 1 
(e~* w ) for the case of an imaginary chemical potential. In particular that implies, for real 
chemical potentials, the impossibility of exploiting the usual even-odd factorization trick 
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Figure 4: Chiral condensate vs. /i 2 at (3 = 1.90. The green (blue) solid line represents the 
polynomial (ratio of polynomials) interpolating function; the dashed lines give the corresponding 
uncertainty, coming from the errors in the parameters of the fit. 

for reducing the number of flavors from 8 to 4. The choice of 8 flavors is therefore linked 
to the need of using an exact Hybrid Monte Carlo algorithm: the last is an unavoidable 
requirement if we want to make a detailed comparison of data at imaginary values of /x 
with data at real values of since systematic effects due to an inexact algorithm could be 
different for the two cases. The choice of a large volume is instead essential if we want to 
make a careful test of the method of analytic continuation, since possible non-analyticities 
will show up only in the thermodynamic limit. The observables we determined are the 
Polyakov loop, the chiral condensate and the fermionic number density (n q ). 

We have considered three values, (3=1.90, 1.45 and 1.30, corresponding to the three 
different regimes exposed in Section and for each we have taken measurements for several 
values of the chemical potential, both imaginary and real. The summary of numerical 
simulations is given in Tables [j], pL 0, 0, [5|, |6|. Simulations have been performed on the 
APE100 and APEmille crates in Bari and on the recently installed computer facilities at 
the INFN apeNEXT Computing Center in Rome. Statistics have been chosen so as to have 
statistical errors well below 1% in most cases: indeed our ability to discern the best among 
a set of possible interpolating functions as well as to detect the exact ranges beyond which 
analytic continuation fails, is strictly related to the statistical precision of our data. 

We have chosen two different strategies for our analysis. We have used the data at 
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Figure 5: As in Fig. [|for the Polyakov loop. 

imaginary chemical potential /U/ to determine the parameters of the interpolating function, 
then we have analytically continued this function to real values of the chemical potential 
and compared there with direct Monte Carlo determinations. In this way we are able to 
test how the method of analytic continuation is able to reproduce the correct physical 
results for real values of fi, and to understand which is the best suited function to do so. 

As an alternative way to analyze our results, we have tried to fit both sets of data 
together, at imaginary and real chemical potential, with several analytic functions and 
using variable ranges for both fi and /i/. In this way, using the x 2 test as a statistical tool, 
we are able to understand in which ranges, if any, the method of analytic continuation makes 
any sense at all, at least within the set of analytic functions taken into considerations. 

In order to fulfill CP invariance, the interpolating function must be a even function 
of jj, for observables, such as the Polyakov loop and the chiral condensate, which do not 
depend explicitly on [i. The fermionic number density, being the logarithmic derivative of 
the partition function with respect to the chemical potential, is instead an odd function of 
ii. 

We separate the discussion of our results for the three different regimes, reflecting the 
different strategies followed in searching for the optimal interpolation and the different 
behaviors observed for the physical observables. 

3.1 The high temperature region [3 > [3e 

For this region we have used two kinds of interpolating functions for the data at imaginary 
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Figure 6: x 2 /d.o.f. of the global fit to real and imaginary chemical potential data for the Polyakov 
loop at [3 = 1.90, for various fitting functions, as the maximum value /Lt max of real chemical potential 
taken into account is varied. Lines connecting data point have been drawn to guide the eye. 



/x: polynomials and ratio of polynomials, the last choice being related to the use of Pade 
approximants suggested in Ref. [45]. For the Polyakov loop and the chiral condensate we 
have considered a second order polynomial in fi 2 , 

A + Bfij + Cfij, (3.1) 

according to the standard approach, and the ratio of two first order polynomials in [i 2 , 

A + Bfi 2 



1 + Cfi 2 



(3.2) 



according to our new proposal. Similarly, for the fermionic number density we have used 
a polynomial of the form 

Afij + Bfi] + Cfi] , (3.3) 

and the ratio 

Our findings at {3 = 1.90 are summarized in Figs. ^, [|, [5| and Table [7|. In Fig. |3] we put 
on the same plot the imaginary part of the fermionic number density as a function of \xi 
and the real part of the fermionic number density as a function of the real /i. The two data 
sets match smoothly at fi = 0, which is a necessary condition for the applicability of the 
method of analytic continuation. The fermionic density approaches two for large values of 
the real chemical potential. This saturation effect is a lattice artifact, which is due to the 
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Figure 7: As in Fig. § for (3 = 1.45. 

fact that no more than two fermions per site can be accommodated on the lattice ("Pauli 
blocking"), and manifests itself for values of the chemical potential which are close to the 
ultraviolet cutoff, i.e. for fx of order 1 (see Refs. [52, 53] for a recent discussion of this 
phenomenon). The solid lines represent the two kinds of interpolating functions, whose 
parameters are determined by a fit on the data at imaginary chemical potential. Here both 
interpolations, polynomial and rational function, nicely reproduce the data at real [i over a 
large interval. Deviations start at values of [i for which the saturation effects are certainly 
important. 

In Fig. |] we show the chiral condensate as a function of \x 2 . Again data at imaginary 
fx, i.e. fi 2 < 0, and data at real /i, i.e. fi 2 > 0, nicely match at fi = 0. This time the 
different behavior of the two kinds of interpolation clearly emerges. The ratio of first order 
polynomials in \i 2 reproduces the data at real /i on a much larger interval than the second 
order polynomial in fj, 2 . Deviations arise for values of real li for which saturation effects 
are probably already important. The same conclusions can be drawn from Fig. [| which 
shows data and interpolations for the Polyakov loop. 

The above conclusions do not change if larger order terms are included in the poly- 
nomial interpolation (|3,1|). In fact, larger order polynomials fail to reproduce the data at 
real ji even earlier in fi than second order polynomials. This is due to the fact that the 
higher order terms of the polynomial are the less accurately determined in the fit to data at 
imaginary fi. On the other side, if in the ratio of polynomials the order of the polynomials 
at the numerator and/or at the denominator is increased, no improvement is observed. 
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Figure 8: As in Fig. for (3 = 1.45. 



It is interesting to compare the parameters in the expansion ( p.3[ ) with those predicted 
by the perturbation theory in /j>/T of the fermionic density. In the infinite temperature 
limit, and taking into account the rotation to real chemical potential in Eq. (3.3), it should 



be A = -nf/N? = -1/2 and B = n f /ir 2 = 8/vr 2 = 0.810569... [13]. These values are in 
rough agreement with our findings at the largest available /3 (see Table @). 

We now expose the results of our combined fits using both sets of data at imaginary 
and real values of fi. The range of values used in the fit is limited on the imaginary 
chemical potential side by the presence of the RW transition, so that we included all data 
with fi] < tt/8. For real chemical potentials, we have considered the possible presence 
of non-analyticities and have repeated our fits for different values of the maximum real 
chemical potential, fi max . 

We report only results obtained for the Polyakov loop: those obtained for the other two 
observables look very similar. We have tried fits with several analytic functions, only a few 
of them being exemplified in Fig. |(| where we report the value of x 2 /d.o.f. as a function of 
Amax- The outcome of our analysis, as evident from Fig. ||, can be summarized as follows: 
acceptable values of x 2 /d.o.f. are obtained once sufficiently higher order polynomials or 
ratio of polynomials are taken into account, but only if /u max is less than about 0.5. Instead 
the value of x 2 /d-o.f. gets sensibly different from one for larger values of // max —* 1) 
regardless of the interpolating function. 

We interpret this result as a proof that, within statistical errors, data at real and 
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imaginary chemical potential can indeed be described by one only analytic function, even 
if in a limited range. Hints of possible non-analyticities appear for real chemical potentials 
fi > 0.5. We believe that the most plausible explanation of them is the onset of saturation 
effects. 

A comment is in order about the use of higher order polynomials. Fig. ^ could give 
the impression that increasing the order of the polynomial sensibly improves the method 
of analytic continuation, since a reasonable x 2 /d.o.i. is obtained for a wider range. This 
could seem to contradict our previous statements about the choice of the optimal function 
for extrapolating data from imaginary values of [a, in fact it is not so. Indeed, one should 
consider that when trying to extrapolate information to fi 2 > having at disposal only 
information from negative values of fi 2 , the use of polynomials in fj, 2 can result in instabil- 
ities in the determination of the coefficients, since a polynomial with positive coefficients 
for /j? < is continued to a polynomial with alternating coefficients for fi 2 > and vice 
versa. These instabilities clearly disappear if data on both sides are available, but of course 
this situation cannot be reproduced for real QCD. 



3.2 The intermediate region fi c < < @e 

Our findings at f3 = 1.45 are summarized in Figs. [7|, ||, ^ and Table |8[ Also in this 
case we have used polynomials and ratio of polynomials as interpolating functions for 
the data at imaginary fi. Here the discussion goes along the same lines as for (3 = 1.90 
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Figure 10: As in Fig. | for /3 = 1.45. 



with one important difference: both in the case of polynomials or ratios of polynomials 
as interpolating functions, and for any of the three observables considered here, a fit at 
imaginary chemical potential with x 2 /d.ol of the order of one is possible only in the interval 
[0, /Zj], with p,j = 0.22 -T- 0.24. This \ii represents the onset of a transition, which shows 
up also as a peak in the chiral condensate susceptibility, centered around that value of jli. 
The fact that the interval in fij available for numerical simulations is shorter makes the 
interpolation, and consequently its continuation, less accurate. Nevertheless, also at this 
(3, the use of ratio of polynomial performs much better than simple polynomials. 

Also in this regime we have performed combined fits using both sets of data at imag- 
inary and real values of [i. In this case the range of values used in the fit is limited, on 
the imaginary chemical potential side, by the presence of the continuation of the physical 
critical line, so that we have included only values fii < 0.20. 

Also in this case we report only results obtained for the Polyakov loop with several 



fitting functions, as exemplified in Fig. 10, as a function of fi ma , x . Results look very similar 
to those obtained for (3 = 1.90, with a few differences: while the quality of the fits obtained 
with the ratio of polynomials does not change with respect to (3 = 1.90, higher order 
polynomials are necessary to obtain reasonable x 2 /d-o.f., and hints of non-analyticities 
show up generally earlier, as a function of /} m ax, than for (3 = 1.90. A possible explanation 
for the different behaviour could reside in the presence of the physical pseudo-critical line 
for imaginary values of fi. 

3.3 The low temperature region (3 < (3 C 

Below f3 c our observables are smooth functions of /}/, with periodicity in fij equal to 7r/4 
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Figure 11: Negative side of the horizontal axis: imaginary part of the fermionic number density 
vs. the imaginary chemical potential at [3 — 1.30. Positive side of the horizontal axis: real part of 
the fermionic number density vs. the real chemical potential at = 1.30. The green (blue) solid line 
represents the Fourier (polynomial) interpolating function and its continuation; the dashed lines 
give the corresponding uncertainty, coming from the errors in the parameters of the fit. 



in the case of the fermionic number density and of the chiral condensate and equal to 
7r/2 in the case of the Polyakov loop. This leads naturally to the use of Fourier sums as 
interpolating functions and, in particular, 



for the chiral condensate, 
for the Polyakov loop and 



A + £cos(8/}/) + Ccos(16/t 7 ) (3.5) 

Acos(4/}i) + £cos(12/}j) (3.6) 
Asin(8/t/) +Ssin(16/t/) (3.7) 



for the fermionic number density, which is odd in ///. 

We summarize our results at (5 = 1.30 in Figs. 11, 12, |13| and Table |9[ The functions 



chosen for our fits, and reported in Eq. ( |3.5| ), ( |3.6j ) and ( |3.7| ), are those containing the 
minimum number of terms necessary to obtain a x 2 /d.o.f. close to one (the use of less 
terms leading to a sensible increase of x 2 /d.o.f.). From Table || it is possible to see that 
the coefficients of the secondary harmonic terms in the Fourier sums (^^) , ([^) and ( |3.7| ) 
are suppressed by a factor of a few tens with respect to the coefficients of the dominant 
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Figure 12: Chiral condensate vs. fi 2 at /3 = 1.30. The green (blue) solid line represents the Fourier 
(polynomial) interpolating function and its continuation; the dashed lines give the corresponding 
uncertainty, coming from the errors in the parameters of the fit. 



harmonic, thus signaling a quite fast convergence of the Fourier sums. It is interesting to 
notice that the term proportional to cos(8/2/) does not appear in Eq. ([O]), i.e. one term in 
the harmonic series for the Polyakov loop seems to be missing. The reason is that in the low 
temperature region center symmetry constrains the Polyakov loop to be zero at fij = ir/8 
(corresponding to the border between the two center sectors), so that all frequencies which 
are even multiples of 4/i/ must be excluded. We have verified that if these frequencies 
are included in the interpolating function, the corresponding coefficients are put to zero 
by the fit. The Fourier sums become sums of hyperbolic sine and cosine functions after 
continuation to real fi, which diverge very rapidly and reproduce only partially the data at 
real (i. The deviation between the extrapolation and the data can be taken as an estimate 
of the pseudo-critical value of This is confirmed by the study of the chiral condensate 
susceptibility, which exhibits a peak centered around that value. 

Alternative attempts with longer Fourier sums or with ratios of Fourier sums did not 
change this scenario. In the case of polynomials as interpolating functions, the behavior 
is similar to Fourier sums if, however, large order polynomials are used (see blue lines in 
Figs. 11 , |l^, 13). In this case, the interpolation of data at imaginary // works in an interval 
shorter than that for Fourier sums. These observations confirm that Fourier sums are 
indeed the natural functions to be used for the analytic continuation in the low temperature 
region. 
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Figure 13: As in Fig. [13 for the Polyakov loop. 



Regarding the global combined fits using both sets of data at imaginary and real values 
of the results of our analysis are reported in Fig. |l4|. In this case we report results for all 
observables, with the same fitting functions, Eq. (|3.5| ), (3.6) and (3/7), used previously. All 



data at imaginary chemical potential are taken into account, since no phase transition at 
all is expected on that side, while data at real \x are limited to a m axi niu m value //-max- 

it 

clearly emerges that, for all observables, both sets of data can be nicely fitted by a common 
analytic function, till /i max reaches the region where the physical pseudo-critical point is 
located; at that point the method of analytic continuation clearly loses any sense. However, 
it is quite interesting to notice that the analytic properties of the partition function for 
imaginary chemical potentials are not influenced at all by the presence of the pseudo-critical 
point at real /i. 

An independent determination of the pseudo-critical chemical potential can be ob- 
tained by the study of the susceptibility of the chiral condensate for real ft, shown in 
Fig. [l^. There is an evident peak at fi ~ 0.28, in good agreement with the determinations 
from the \ 2 test method. 



4. Conclusions and outlook 



We have studied the method of analytic continuation in a theory which does not suffer from 
the sign problem and have looked for better interpolating functions at imaginary fj,, to be 
used instead of the polynomial, as has been done in most cases so far in the literature. 
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Figure 14: x 2 /d.o.f. of the global fit to real and imaginary chemical potential data for the various 
observable at = 1.30, as a function of the maximum value /x max of real chemical potential taken 
into account. Fitting functions are the same used for Figs, [ll], [lj, |l3| and reported in Eq. (3.5), 
(3.6) and (3/7). The values of x 2 /d.o.f. are quite stable around one for all observables, till /2 max 
crosses the pseudo-critical point. 



We have verified that data at real and imaginary chemical potential can indeed be 
well described by common suitable analytic functions, in appropriate ranges, and we have 
found that a considerable improvement can be achieved, when extrapolating data from 
imaginary to real chemical potentials, if ratios of polynomials (or equivalently Pade ap- 
proximants [45]), are used as interpolating functions, if the temperature is larger than the 
pseudo-critical one at zero chemical potential. Below that value, instead, Fourier sums 
seem to be the best Ansatz, as expected and tested also in other contexts [33,34]. 

The deviations from analyticity and between the extrapolated functions and the data 
at real chemical potential have different explanations, according to the temperature regime. 
Above the temperature of the RW endpoint they arise most likely from unphysical satu- 
ration effects, due to the lattice discretization ("Pauli blocking"). In the intermediate 
regime, deviations stem also from the limited range of the interval in the imaginary chem- 
ical potential for the numerical simulations, which makes the interpolation less easy: this 
is caused by the presence of a pseudo-critical point for imaginary values of the chemical 
potential, which could also contribute to restrict the range where analytic continuation 
can be applied. Finally, in the low temperature regime, deviations necessarily appear in 
correspondence of the transition at real chemical potential. 

The lessons we have learned from this study and which could be applied to the physi- 
cally interesting case of SU(3) can be summarized as follows: 

• above the pseudo-critical temperature, ratio of polynomials should be used as interpo- 
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Figure 15: Susceptibility of the chiral condensate at f3 = 1.30 for real ft. 



lating functions instead of polynomials; their continuation to real chemical potentials 
is the more reliable the larger is the interval of imaginary chemical potential where 
they succeed in interpolating data; 

• below the pseudo-critical temperature, one should surely use Fourier sums: they 
nicely reproduce data at imaginary chemical potentials, but are extrapolated to hy- 
perbolic functions which rapidly diverge at real chemical potentials; nevertheless, 
analytic continuation works fairly well till the pseudo-critical value of the real chem- 
ical potential is reached. 
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Table 1: Summary of the simulations at j3 — 1.90 and real chemical potential, (i = £lr. 
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Table 2: Summary of the simulations at {3 = 1.90 and imaginary chemical potential, fi = ifu. 
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Table 3: Summary of the simulations at (3 = 1.45 and real chemical potential, fi = fin. 
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Table 4: Summary of the simulations at (3 = 1.45 and imaginary chemical potential, /t = i/t/. 
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Table 5: Summary of the simulations at j3 = 1.30 and real chemical potential, fi = flu. 
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Table 6: Summary of the simulations at j3 = 1.30 and imaginary chemical potential, (l = 
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Table 7: Parameters of the interpolations of imaginary chemical potential data at (3 = 1.90. 
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0.25072(18) 


-1.200(23) 


-10.76(43) 


0.79 


0.24 


Q{n q ) 


Eq. 


M> 


-0.36813(79) 


3.716(42) 


-2.20(27) 


1.33 
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Table 8: Parameters of the interpolations of imaginary chemical potential data at (3 = 1.45. The 
last column gives the largest value of /2/ included in the fit. 



observable 


function 


A 


B 


C 


X 2 /d.o.f. 




Eq. (p|) 
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1.25 


{•H>) 


Eq. (p|) 
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Eq. (pD 


0.12426(8) 


0.002238(81) 




1.04 



Table 9: Parameters of the interpolations of imaginary chemical potential data at (3 — 1.30. 



